function[x] = quadsolver(p, q)
	d = (p/2)^2 - q;
	if(d < 0)
		x = [];
	else
		x(1) = -p/2 - sqrt(d);
		x(2) = -p/2 + sqrt(d);
	endif
endfunction